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Electrochemical control of quantum interference in anthraquinone-based molecular 

switches 
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Using first-principles calculations we analyze the electronic transport properties of a recently 
proposed anthraquinone based electrochemical switch. Robust conductance on/off ratios of sev- 
eral orders of magnitude are observed due to destructive quantum interference present in the an- 
thraquinone, but absent in the hydroquinone molecular bridge. A simple explanation of the inter- 
ference effect is achieved by transforming the frontier molecular orbitals into localized molecular 
orbitals thereby obtaining a minimal tight-binding model describing the transport in the relevant 
energy range in terms of hopping via the localized orbitals. The topology of the tight-binding model, 
which is dictated by the symmetries of the molecular orbitals, determines the amount of quantum 
interference. 

PACS numbers: 



I. INTRODUCTION 

Molecular electronics is a very active area of research 
which holds the promise of continuing the miniaturiza- 
tion of active devices beyond the limits of standard silicon 
technologies [ll-Q- Of particular interest are molecules 
which can be switched between two distinct states with 
different electronic conductance Q • Reversible switching 
between the 'on' and 'off' states can be mediated with 
light 0^31: bias voltage M, or by changes in the elec- 
trochemical environment \\SX \Va . It has recently been 
proposed by van Dijk et al. [12| to use an anthraquinone 
based molecule as an electrochemical switch. Cyclic 
voltamctry showed that the anthraquinone (AQ) could 
be revcrsibly reduced to hydroquinone (HQ), with a sig- 
nificantly different UV-vis absorption spectra. In this 
communication we calculate the conductance properties 
of the proposed switch and show that very large on/off 
ratios should be expected due to destructive quantum 
interference effects in the off-state (AQ) molecule. 

It has previously been shown that quantum interfer- 
ence effects can have a dramatic impact on the electronic 
transport through a single molecule |13l - [la |. As an ex- 
ample, experiments have shown that benzene has very 
different transport properties depending on the position 
of the anchor groups IITI llSj . This has been confirmed 
by theoretical works [19[ and explained in terms of in- 
terfering pathways. Interference effects in benzene have 
also been studied for the non-coherent coulomb block- 
ade regime [20, [2l[ . Related theoretical studies has re- 
cently been concerned with interference effects in ni- 
trobenzene 22| and in cross conjugated molecules |23l.[24| 
showing transmission anti-resonances close to the Fermi 
energy, thus potentially making the destructive quantum 
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interference features useful in molecular electronics ap- 
plications. Contrary to the previously discussed cross 
conjugated molecules, the AQ-HQ setup offers the pos- 
sibility to electrochemically switch between a situation 
with and without transmission anti-resonances. 

We show in this work that the interference effects in an- 
thraquinone electrochemical switches can be understood 
from a minimal tight-binding (TB) model constructed 
from the frontier molecular orbitals of the free molecule. 
Such simple models are clearly desirable for an intuitive 
understanding and for a fast screening of molecules ex- 
hibiting interesting interference features. 



II. DFT-LANDAUER TRANSMISSIONS 

We calculate the electronic conductance of the molec- 
ular contacts shown in Figure [T] (a)-(b), using a combi- 
nation of density functional theory (DFT) and Green's 
function (GF) methods. In Figure [T] (a) the molecule 
is terminated with two S-atoms that bind to both Au 
surface. This contact resembles the situation in e.g. me- 
chanical break junction experiments. Figure [T](b) repre- 
sents an STM experimental situation, where the molecule 
is only bound to one Au surface, while the other contact 
is an STM tip. 

For both setups, we initially relax the molecule and the 
two closest Au layers are using a double-zeta with polar- 
ization (dzp) atomic basis set [25| using the grid-based 
projector augmented wave method GPAW [2g] with 
the Perdew-Burke-Ernzerhof GGA exchange-correlation 
functional ^\ . In the relaxed configuration, the S-atoms 
bind to Au at a bridge site slightly shifted toward the 
hollow site. Following the standard DFT-Landauer ap- 
proach as described in Ref. [2^, we calculate the zero- 
bias transmission function, T{E). The low-bias conduc- 
tance can finally be obtained from the Landauer formula, 
G ~ (2e^//i) T{Ep), where Ep is the Fermi energy. 
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Figure 1; Atomic structure of HQ molecule coupled to gold 
electrodes (a). The two hydrogens in the circles can reversibly 
be attached or removed electrochemically forming HQ (with 
H-atoms) and AQ (without H-atoms). Panel (b) shows an 
STM setup where the molecule is bound to only one An sur- 
face. Panels (c) and (d) show the electronic transmission func- 
tions for the two molecular contacts in panel (a) and (b), re- 
spectively. Solid curves correspond to HQ while dashed curves 
show AQ results. Both setup give large on-off ratios at the 
Fermi level. 



Figure [T](c) shows the electronic transmission function 
for the AQ (dashed black) and HQ (solid red) molecu- 
lar bridges show in panel (a). There is a strikingly large 
difference of six orders of magnitude in the transmission 
around the Fermi energy. While the HQ has a relatively 
large transmission of around 0.01, the AQ displays a clear 
transmission anti-resonance with a minimum transmis- 
sion of 10"^. 

Figure [T] (d) shows the transmission functions of the 
HQ and AQ in the STM setup (Figure [2 (b)). The 
transmission values are generally lower by a factor 100 
compared to the contact in Figure [T] (a). This is due 
to the weak tunneling contact between the molecule and 
the tip. We again observe a clear transmission dip for the 
AQ but not for the HQ. Both the transmission peaks and 
dip are shifted toward lower energies in the STM setup, 
and the antiresonance is no longer right at the Fermi 
level. We attribute this to the different coupling to the 
Au contact which cause different shifts of the molecular 
levels. It might still be possible that even the trans- 
mission antiresonance away from the Fermi level can be 
observed experimentally by tuning the molecular levels 
with the electrochemical gate voltage. The qualitative 
agreement between the two different contact situations 
show that the interference feature in the AQ is an intrin- 
sic feature of the molecule, which is only quantitatively 



affected by the coupling to the metal. The large on/off 
ratio should therefore be a robust feature of the AQ-HQ 
switch. In a recent STM conductance measurement on 
quinone-oligo(phenylene vinylene) (Q-OPV) molecules in 
an electrochemically controlled environment [29|, it was 
found that the reduced HQ-OPV molecule has a 40 times 
larger conductance than the Q-OPV. This results is in 
qualitative agreement with our finding that HQ has a 
larger conductance than AQ. 

It is well known that DFT in general does not pro- 
vide an accurate description of molecular energy gaps 
and level alignment at metallic surfaces [30, |31|. In order 
to verify that the interference effect is not sensitive to 
changes in the molecular level positions we have changed 
"by hand" the molecular orbitals using a scissors operator 
technique similar to that described in Ref. [S^l ■ Shifting 
the energies of the occupied states down by 1 cV and 
the energies of the unoccpicd states up by 1 eV slightly 
shifts the transmission minimum of the AQ, but does not 
change the overall picture. The large on/off ratio should 
thus be a rather robust feature of the AQ-HQ switch. 

In the remainder of the paper we analyze the AQ-HQ 
switch in terms of two simplifying models that offers com- 
plementary insights to the underlying physics of the in- 
terference effects. 



III. HUCKEL TT-MODEL 

In order to gain insight in the physics underlying the 
interference effects, we describe the AQ/HQ molecule 
with a Hiickel 7r-electron model, in which there is one 
Pz-orbital at every carbon atom. We assume only near- 
est neighbour interaction with a hopping parameter, 
t = — 3eV, and constant on-site energy Eq = OeV. In 
the oxidized AQ state, we further assume that the oxy- 
gens contribute with one Pz orbital each with the same 
on-site and hopping parameters as the C atoms. In the 
reduced HQ state this extra electron will be passivated by 
the hydrogens, and does not contribute to the tt system. 
The schematics of the central parts of the AQ and HQ 
molecules in the tt- model arc shown as insets in Fig. [2] (a) 
and (b), respectively. The transmission in the 7r-model 
is calculated using a wide band limit approximation for 
the electrode self-energies, which adds an imaginary part, 
ij, to the on-site energy of the C-atoms connected to the 
sulfurs. The sulfur atoms themself are not included in 
the model. 

Figure [2] shows the full DFT transmission (solid black) 
together with the 7r-modcl transmission (dashed red) for 
the AQ (a) and HQ (b) molecules. Both the transmission 
anti-resonance in the AQ and the relatively large trans- 
mission for the HQ found in the DFT-GF calculations 
are well reproduced by the 7r-model. The good agree- 
ment between the DFT-GF and the 7r-model results is 
not obvious as it might be expected that both the on- 
site energy of and the hopping to the oxygens are quite 
different from the values for pure carbon. However, it is 
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Figure 2: Transmission of the AQ (a) and HQ (b) molecules 
calculated with a Hiickel 7r-model (dashed red) and compared 
with the DFT-GF results (solid black). The insets show the 
structure used in the Hiickel model 



states. Often a qualitatively correct description can be 
obtained by considering only the single molecular level 
closest to Fermi energy, i.e. either the highest occu- 
pied molecular orbital (HOMO) or the lowest unoccupied 
molecular orbital (LUMO). Assuming a constant density 
of states (DOS) in the electrodes (the so-called wide band 
limit), the transmission function of a single- level model 
becomes a Lorentzian centered at the energy of the MO 
and with a width, F determined by the electrode DOS 
and the spatial overlap of the MO with the states in the 
electrode. The frontier MOs of the free (i.e. not coupled 
to gold) AQ and HQ are shown in Figure [31 together with 
the level energies. 
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encouraging that the simple 7r-model reproduce the DFT 
results so well. The almost exact coexistence of the trans- 
mission minimum at the Fermi energy, might not be a 
general result, since the exact position of the DFT trans- 
mission minimum is sensitive to the bonding geometries, 
as illustrated by the differences in transmission between 
the two setups in Fig. [TJ 

The results of the 7r-model shows that switching from 
HQ to AQ can be thought of as effectively adding an or- 
bital to the TT-system. An electron in AQ travelling from 
left to right can either do it directly or via the extra side 
orbital, with the two pathways intefering destructively at 
the Fermi energy. In the HQ, on the other hand, there 
is only the direct path and no interference effects should 
be expected. 



IV. LOCALIZED MOLECULAR ORBITALS 

While the Hiickel model succesfully explains the in- 
terference effects in terms of an extra 7r-orbital in the 
oxidized AQ state, it is desirable to understand the in- 
terference effects in terms of the frontier orbitals. To 
this end, we develop in this section a simple and general 
method for analyzing quantum interference in molecu- 
lar junctions using the AQ/HQ switch in Figure [1] (a) 
as a specific example. In the physically relevant energy 
range around the Fermi level, the transmission through 
a molecular junction is governed by the frontier molecu- 
lar orbitals (MOs) and their overlap with the electrode 



Figure 3: Molecular eigenstates of the isolated AQ and HQ, 
respectively. Note that the HOMO and LUMO for the HQ 
have similar shapes as the LUMO and LUMO-l-1 for the AQ. 
The effect of the two extra H atoms in the HQ is thus merrily 
to occupy one extra orbital without significantly changing the 
orbital shapes. The energies of the eigenstates relative to the 
HOMO levels are given in units of eV. 

In Figure [l] (c) the transmission function of HQ can 
approximately be described by a sum of two Lorentzians 
centered at the HOMO and LUMO energies, respectively. 
In contrast this is not the case for AQ. The pronounced 
difference in the transmission properties of AQ and HQ 
is difficult to detect from the shape and energies of the 
frontier MOs (see Figure ^. The difference, however, 
becomes clear if we transform the frontier MOs, ipm, into 
localized molecular orbitals (LMOs), 0„: 



^n / ^ ^ nin'Yn 



(1) 



The expansion coefficients are determined following the 
standard procedure for constructing maximally localized 
Wannicr functions (for solids) or maximally localized 
molecular orbitals (for molecules) [33, |3J| . We thus seek 
to minimize the spread of the LMOs measured by the 
sum of their second moments: 
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Figure 4: Localized molecular orbitals (LMOs) and the corresponding TB models. The LMOs have been constructed from (a) 
the HOMO-l,HOMO, and LUMO of AQ (b) the HOMO and LUMO of AQ, and (c) the HOMO and LUMO of HQ. TB models 
(a) and (b) offers the opportunity for destructive interference and transmission anti-resonances due the multiple paths through 
the molecule. 



The mapping from the molecular eigcnstates (■0m) to 
the LMOs {4>n) described by the unitary matrix. U. is 
uniquely determined by the requirement that the result- 
ing LMOs should be maximally localized as defined by 
Eq. ((21). Usually, one maps the entire set of occupied 
MOs into an equivalent set of LMOs. In the present case, 
however, we are interested in describing the physics in the 
vicinity of the Fermi level and thus we transform only the 
MOs closest to Ep. Using the LMOs as basis states we 
can construct an effective TB-like Hamiltonian for the 
molecule. As an example Figure |4|a) shows the LMOs 
obtained from the HOMO-1, HOMO, and LUMO of the 
AQ. The corresponding TB Hamiltonian is obtained as 



H 



U^diag(eH-i,£H,£L)U 



Pl £i Pb^ 



(3) 



where eH-i,£H,£L denote the molecular energy levels. 
In the considered case the TB Hamiltonian has certain 
symmetries which derives from symmetries of the LMOs: 
The two bottom LMOs of Figure IHa) both have on- 
site energy Eq = 0.07 eV (relative to the AQ HOMO 
level) but are localized on the left and right parts of the 
molecule, respectively. The topmost LMO has on-site 
energy e\ = 1.56 eV and is locahzed at the center of the 
molecule with a vanishing weight at the terminating S 
atoms. Assuming a wide band approximation, the lead 
self-energies only adds an imaginary part to the left and 
right on-site energies, Eq ~^ £o + il- This express the fact 
that the molecule is only coupled to the leads via the two 
LMOs located at the ends of the molecule. The topology 
of the effective three-site TB model is sketched in the 
upper panel of Figure S] (a). The hopping parameters 
are a = 0.07 eV and Pl = Pr = 0.25 eV. It is intuitively 
clear that the two pathways through the molecule, i.e. 
hopping from the left to the right LMO either directly or 
via the central LMO, can lead to interference effects in 



the transmission. 

The LMOs and corresponding TB model obtained 
through the procedure described above depends on the 
initial set of MOs used for the mapping. In Figure |H[b) 
we show the LMOs and the corresponding TB topol- 
ogy obtained when the mapping is performed from the 
HOMO and LUMO of the AQ. The obtained TB param- 
eters are eq = 0.15 eV, ei = 1.56 eV, /3 = 0.35 eV. As for 
the three-state case one can anticipate interference effects 
due to transmission though two different pathways - ei- 
ther direct hopping or hopping via the LMO localized at 
the center. 

Figure [5] (a) shows the transmission of the AQ (on a 
log scale) calculated using the full DFT-GF setup (dot- 
ted black), which is the same as shown in Figure [T] 
(c). Also shown are the transmissions obtained from 
the three- and two-site TB models in Figure S] (a)-(b). 
As we do not have any reference to the electrode Fermi 
energy in the TB models, we have aligned the LUMO 
transmission peak (at E = 0.6 eV) of the TB transmis- 
sions to the DFT-GF result. Evidently, the three-site 
model (solid red) qualitatively reproduces the transmis- 
sion anti-resonance and the overall shape of the trans- 
mission function is better reproduced by the three-site 
model than the two-site model (dashed blue), which in 
general yields a too large transmission and the transmis- 
sion anti-resonance is too high in energy. It can readily 
be found that the two-site model has a transmission zero 
at energy E = ii, whereas the three-site model has the 
anti-resonance at £■ = ei — PLPRJa. Note that in both 
models the anti-resonance energy is independent of the 
on-site energies, eo and eoj and of the LMOs couphng 
to the electrodes. It is not surprising that the three-site 
model gives a better description of the full calculation 
as the energies and wavefunctions of the HOMO and 
HOMO-1 are very similar, see Figure [31 We note that 
the quantitative differences between the full and model 



calculations are due to (i) neglect of all but a few MOs 
in the construction of the LMOs, and (ii) the wide band 
approximation. 
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Figure 5: (a) Transmission function of the AQ coupled to 
gold electrodes calculated using the DFT-GF method (dot- 
ted black) and for the two-site (dashed blue) and three-site 
(solid red) TB models in Figure [1(b) and (c). (b): DFT- 
GF transmission (dotted black) for the HQ together with TB 
transmission for the two-site model (dashed blue) in Figure 
[4] (c) and for a four site-model including two extra occupied 
MOs (solid red). 

We now turn to the HQ molecule. Figure ID^c) shows 
the LMOs constructed from the HOMO and LUMO or- 
bitals, and the corresponding TB model. The on-site 
energies are Eq = 0.55 eV, £i = 0.97 eV (relative to the 
HQ HOMO level), and fi = -0.73 eV. The topology of 
the HQ two level model is qualitatively different from the 
AQ topologies in that there is only one possible pathway 
through the molecule and thus interference effects are 
not anticipated. Figure [5{b) compares the transmission 
through the HQ obtained from the full calculation (dot- 
ted black), the two-site TB model in Fig. \M,c) (dashed 
blue), and a four-site model including also the HOMO-2 
and HOMO-1 (solid red). Both TB models reproduce 
the full result, and neither of them shows interference 
around the Fermi level. The four-site model gives an 
anti-resonance at energies below the HOMO level which 
correspond to the anti-resonance in the AQ. 



A. Application of LMOs to a different molecule 

To further validate the simple model calculations us- 
ing LMOs of the isolated molecules, we consider the 
cross-conjugated molecule shown in Figure [6] (b). The 
transmission properties of similar molecules have pre- 
viously been analyzed [35| and clear transmission anti- 



resonances were observed, similar to the AQ transmis- 
sion shown above. In Figure |6] (a) we compare the 
transmission through the cross-conjugated molecule at- 
tached to Au(lll) electrodes as above and calculated 
with DFT-GF (dotted black), and with two TB mod- 
els. Our calculated transmission agrees with the results 
in Ref . [Sa] ■ The three-site model (solid red) gives again a 
qualitatively correct description of the transmission anti- 
resonance. The corresponding LMOs, constructed from 
the H0M0-1,H0M0, and LUMO, are shown in Figure [6] 
(e)-(g). The two-site orbitals are shown in Figure |6] (c)- 
(d). As for the AQ, the two-site model yields a transmis- 
sion anti-resonance higher in energy than both the three- 
site model and the full calculation. The good agreement 
between the three-site TB model and the fuU DFT-GF 
calculation further illustrates that the simple minimal 
TB model provides a transparent and qualitative correct 
description of the full transport problem. The molecu- 
lar level energies for the free cross-conjugated molecules 
are Ehomo-i = -0.4 eV, enoMo = 0.0 eV, Elumo = 2.0 eV, 
relative to the HOMO level. Contrary to the AQ, the 
HOMO-1 and HOMO levels are clearly separated in en- 
ergy, and the better results obtained with the three- 
site model including the HOMO-1 orbital indicates that 
quantum interference effects in molecules more generally 
involves the HOMO-1, HOMO, and LUMO orbitals and 
not only the HOMO and LUMO, as might be anticipated. 
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Figure 6: (a): Transmission through the molecule shown 
in panel (b) calculated with the DFT-GF method with the 
molecule coupled to gold electrodes (solid red) and within 
the wideband approximation using the the two- (dashed blue) 
and three-site (dotted black) models in Figure|4](b) - (c). The 
corresponding LMOs are shown in panel (c)-(d) and (e)-(g) 
for the two- and three site models, respectively. 



V. SUMMARY 

In summary, we have perform.ed ab initio transport 
calculations for a recently proposed anthraquinone based 
electrochemical molecular switch [i2|- Due to destruc- 
tive quantum interference effects the conductance in the 
reduced ('off') state is dramatically suppressed by sev- 
eral orders of magnitude as compared to the oxidized 
('on') state. The desctructive interference were explained 
within a simple Hiickel 7r-model by the addition of two 
extra pz orbitals on the oxygen atoms in the oxidized 
anthraquinone state. The interference effects were fur- 
ther rationalized by mapping the full transport problem 
to an effective minimal TB model obtained via a trans- 



formation of the frontier molecular orbitals into localized 
molecular orbitals. The approach can be used to cat- 
egorize any given molecular junction according to the 
topology of the effective TB model and offers a trans- 
parent and fast method for computational screening for 
molecules exhibiting interference features. 
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